Spherical Scalar Field Halo in Galaxies 



Tonatiuh Matos, 1 F. Siddhartha Guzman 1 * and Dario Nunez 2 
1 Departamento de Fisica, Centro de Investigation y de Estudios Avanzados del IPN, A. P. 14--740, 07000 Mexico D.F., 

MEXICO. 

2 Instituto de Ciencias Nucleares, Universidad National Autonoma de Mexico, A. P. 70-543, 04510 Mexico D.F., MEXICO 

(February 1, 2008) 

We study a spherically symmetric fluctuation of scalar dark matter in the cosmos and show that 
it could be the dark matter in galaxies, provided that the scalar field has an exponential potential 
whose overall sign is negative and whose exponent is constrained observationally by the rotation 
velocities of galaxies. The local space-time of the fluctuation contains a three dimensional space-like 
hypersurface with surplus of angle. 

PACS numbers: 95.35.+d, 95.35.G. 

The existence of dark matter in the Universe has been firmly established by astronomical observations at very dif- 
ferent length-scales, ranging from single galaxies, to clusters of galaxies, up to cosmological scale (see for example 0). 
A large fraction of the mass needed to produce the observed dynamical effects in all these very different systems is not 
seen. At the galactic scale, the problem is clearly posed: The measurements of rotation curves (tangential velocities 
of objects) in spiral galaxies show that the coplanar orbital motion of gas in the outer parts of these galaxies keeps a 
more or less constant velocity up to several luminous radii Q , forming a radii independent curve in the outer parts of 
the rotational curves profile; a motion which does not correspond to the one due to the observed matter distribution, 
hence there must be present some type of dark matter causing the observed motion. The flat profile of the rotational 
curves is maybe the main feature observed in many galaxies. It is believed that the dark matter in galaxies has an 
almost spherical distribution which decays like 1/r 2 . With this distribution of some kind of matter it is possible to fit 
the rotational curves of galaxies quite well Q. Nevertheless, the main question of the dark matter problem remains; 
which is the nature of the dark matter in galaxies? The problem is not easy to solve, it is not sufficient to find out 
an exotic particle which could exist in galaxies in the low energy regime of some theory. It is necessary to show as 
well, that this particle (baryonic or exotic) distributes in a very similar manner in all these galaxies, and finally, to 
give some reason for its existence in galaxies. 

In previous works it has been explored, with considerable success, the possibility that scalar fields could be the dark 
matter in spiral galaxies by assuming that the scalar dark matter distributes as an axially symmetric halo ^J^] . The 
idea of these works is to explore whether a scalar field can fluctuate along the history of the Universe and thus forming 
concentrations of scalar field density. If, for example, the scalar field evolves with a scalar field potential U(<£>) ~ 3> 2 , 
the evolution of this scalar field will be similar to the evolution of a perfect fluid with equation of state p = 0, i.e., it 
would evolve as cold dark matter ||. However, it is not clear whether a spherical scalar field fluctuation can serve as 
dark matter in galaxies. In this letter we show that this could be the case. We assume that the halo of a galaxy is 
a spherical fluctuation of cosmological scalar dark matter and study the consequences for the space-time background 
at this scale, in order to restrict the state equation corresponding to the dark matter inside the fluctuation. We start 
from the general spherically symmetric line element and find out the conditions on the metric in order that the test 
particles in the galaxy possess a flat rotation curve in the region where the scalar field (the dark matter) dominates. 
Finally we show that a spherical fluctuation of the scalar field could be the dark matter in galaxies. 

Assuming thus that the dark matter is scalar, we start with the energy momentum tensor T^ u = <& tli <& tU — 
l/2g AI „<F <J <I> iCr — g M „U($), being $ the scalar field and V(<&) the scalar potential. The Klein-Gordon and Einstein 
equations respectively are: 

**-^ = 

where i? M „ is the Ricci tensor, \J—g the determinant of the metric, kq = 8irG and a semicolon stands for covariant 
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derivative according to the background space-time; /x, v — 0, 1, 2, 3. 



Assuming that the halo has spherical symmetry and that dragging effects on stars and dust are inappreciable, i.e. 
the space-time is static, the following line element is the appropriate 

ds 2 = -B(r)dt 2 + A{r)dr 2 + r 2 d6 2 + r 2 sin 2 8d<p 2 (1) 

where A and B are arbitrary functions of the coordinate r. Following the analysis made for axisymmetric stationary 
space-times J?J , we consider the Lagrangian for a test particle travelling on the space time described by ([!]) which is 

2C = -Bi 2 + Ar 2 + r 2 6 2 + r 2 sin 2 8tp 2 (2) 

where a dot means derivative with respect to the proper time. From (^|) the generalized momenta read: 

p t = -E= -Bi (3) 

Pr = Ar (4) 

Pe = L e = r 2 6 (5) 

Pip = L v = r 2 sin 2 Qtp (6) 

being E the total energy of a test particle and Li the component of its angular momentum. It can be defined the 
Hamiltonian 7i — p^q^ — C and after rescalling the proper time for the lagrangian to equal 1/2 for time-like geodesies, 
the geodesic equation for material particles (stars and dust) arises 

L 2 

being L\ — L 2 , + iHhj the first integral corresponding to the squared total angular momentum. We are interested in 
circular and stable motion of test particles, therefore the following conditions must be satisfied 

i) r = 0, circular trajectories 

ii) 9 ^ — 0, extreme ones 

\extr > 0, and stable. 

being V(r) = [l + L^/r 2 — E 2 /B] /A. Following || it is found that the tangential velocity of the test particle is 



rB' 

tangential _ (p _ fg\ 

V 2B y ' 

where ' means derivative with respect to r. It is easy to show that if flat rotation curves are required, it arises the 
following flat curve condition from (||), that is B = B$r l with I = 2(v v ) 2 . With the flat curve condition, metric (|l|) 
becomes 

ds 2 = -B r l dt 2 + A(r)dr 2 + r 2 d0 2 + r 2 sin 2 6d(p 2 (9) 

This result is not surprising. Remember that the Newtonian potential tp is defined as goo = —exp(2ip) = —X — 2ip . 

On the other side, the observed rotational curve profile in the dark matter dominated region is such that the rotational 
velocity v v of the stars is constant, the force is then given by F — —(v v ) 2 /r, which respective Newtonian potential is 
i> = (v v ) 2 ln(r). If we now read the Newtonian potential from the metric (||), we just obtain the same result. Metric 
) is then the metric of the general relativistic version of a matter distribution, which test particles move in constant 
rotational curves. Function A will be determined by the kind of substance we are supposing the dark matter is made 
of. Assuming the flat curve condition in the scalar dark matter hypothesis, we are in the position to write down the 
set of field equations. Using (0), the Klein Gordon equation reads 
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and the Einstein equations are 
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(11) 
(12) 

(13) 



In order to solve equations (|Tl l3), observe that the combination of the previous equations [(2 — i)(eq. (O)— 4 (eq. 
H)+(2 + 0(eq. Pb] implies 
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«o(2 - r 2 



(14) 



This is a very important result, namely the scalar potential goes always as 1/r 2 for a spherically symmetric metric 
with the flat curve condition. It is remarkable that this behavior of the stress tensor coincides with the expected 
behavior of the energy density of the dark matter in a galaxy. We can go further and solve the field equations, the 
general solution of equations ( |IT[jT^ ) is 

A{r) = (4 - I 2 ) I (4 + C (4 - I 2 ) r-('+ 2 )) , 

being C an integration constant and we can thus integrate the function Nevertheless, in this letter we consider 
the most simple solution of the field equations with C = 0. Observe that for this particular solution the stress tensor 
goes like 1/r 2 . The energy momentum tensor is made essentially of two parts. One is the scalar potential and the 
other one contains products of the derivatives of the scalar field, both going as 1/r 2 . Furthermore, as (Q, r ) 2 ~ l/?' 2 , 
this means that <E> ~ ln(r), implying that the scalar potential is exponential V ~ exp(2a$) such as has been found 
useful for structure formation scenarios ^|Jl0|| and scaling solutions with a primordial scalar field in the cosmological 
context P-^2[ including quintessential scenarios Jl3|]. Thus, the particular solution for the system (|l^ - |i~3| ) that we 
are considering is 
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exp[-2W^($-$ )] 



(15) 
(16) 
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where ([[5]) and ( |l6| ) approach asymptotically (r — > 00) the case with m = 2, n — 21/(2—1) in the general study 
of the global properties of spherically symmetric solutions in dimensionally reduced space-times ]T^| . Function A 
corresponds to an exact solution of the Einstein equations of a spherically symmetric space-time, in which the matter 
contents is a scalar field with an exponential potential. Let us perform the rescalling r 2 — > 4r 2 / (4 — Z 2 ). In this case 
the three dimensional space corresponds to a surplus of angle (analogous to the deficit of angle) one; the metric reads 



ds 2 



-B Q r l dt 2 + dr 2 + 



A -I 2 



[d9 2 



(18) 



for which the two dimensional hypersurface area is 47rr 2 x 4/ (A — I 2 ) = Anr 2 /(l — (v v ) A ). Observe that if the rotational 
velocity of the test particles were the speed of light v v — * 1, this area would grow very fast. Nevertheless, for a typical 
galaxy, the rotational velocities are v v ~ 10~ 3 (300fcm/s), in this case the rate of the difference of this hypersurface 
area and a flat one is (v v ) A /(l — {v v ) A ) ~ 10~ 12 , which is too small to be measured, but sufficient to give the right 
behavior of the motion of stars in a galaxy. 

Let us consider the components of the scalar field as those of a perfect fluid, it is found that the components of the 
stress-energy tensor have the following form 
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I 2 I 

P = T V = fti> ' (20) 

while the angular pressures are Pg = P v = —p. The analysis of an axially symmetric perfect fluid in general is given 
in where a similar result was found (see also [0). 

The effective density ( |l9| ) depends on the velocities of the stars in the galaxy, — p = (u v ) 4 /(l — {v v ) 4 ) x l/(re r 2 ) 
which for the typical velocities in a galaxy is — p ~ 10~ 12 x l/(n^r 2 ) 1 while the effective radial pressure is 
P = (v v ) 2 ((v ip )' 2 + 2)/(l - (v v ) 2 ) x l/(K Q r 2 ) ~ 10~ 6 x l/(K r 2 ), i.e., six orders of magnitude greater than the 
scalar field density. This is the reason why it is not possible to understand a galaxy with Newtonian dynamics. 
Newton theory is the limit of the Einstein theory for weak fields, small velocities but also for small pressures (in com- 
parison with densities). A galaxy fulfills the first two conditions, but it has pressures six orders of magnitude bigger 
than the dark matter density, which is the dominating density in a galaxy. This effective pressure is the responsible 
for the behavior of the flat rotation curves in the dark matter dominated part of the galaxies. 



Metric (|18|) is not asymptotically flat, it could not be so. An asymptotically flat metric behaves necessarily like a 
Newtonian potential provoking that the velocity profile somewhere decays, which is not the observed case in galaxies. 
Nevertheless, the energy density in the halo of the galaxy decays as 

-p - = -^-Porit (21) 

where Hq 1 — VS/h 10 6 Kpc is the Hubble parameter and p cr it is the critical density of the Universe. This means that 
after a relative small distance r cr n ~ y/3/h 2 « 3Kpc the effective density of the halo is similar as the critical density 
of the Universe. One expects, of course, that the matter density around a galaxy is smaller than the critical density 
p5[ , say paround ~ 0-06p cr it, then r cr it ~ \AKpc. Observe also that metric ( |l8| ) has an almost flat three dimensional 
space-like hypersurface. The difference between a flat three dimensional hypersurface area and the three dimensional 
hypersurface area of metric (Q) is ~ 10 -12 , this is the reason why the space-time of a galaxies seems to be so flat. 
We think that these results show that it is possible that the scalar field could be the missing matter (the dark matter) 
of galaxies and maybe of the Universe. 

Possibly the greatest problem with the present model is the physical origin of the exponential potential ([l7]) . Firstly, 
its sign is necessarily opposite to that of the exponential potentials that have been considered in quintessence cosmolo- 
gies |9 - 11 13fl. Secondly, although exponential scalar potentials with an overall negative sign do arise from dimensional 



reduction of higher-dimensional gravity with the extra dimensions forming a compact Einstein space of dimension 
n > 2 such models also constrain the parameter I = 2n/ (n + 2) to take values I > 1, which are inconsistent with 
its interpretation as I — 2 (u^) 2 for velocities v v of the order of magnitude of the rotation velocity of galaxies. Similar 
considerations apply to the exponent of a single exponential potential obtained by the dimensional reduction of a 
theory with a higher-dimensional cosmological constant and Ricci-flat internal space [^6|. Nonetheless, exponential 
potentials can arise in a variety of ways in stringy gravity, possibly via symmetry breaking or similar mechanisms, 
and so we are hopeful that a natural origin can be found for potentials of the type considered here. 
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